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We consider disk accretion resulting purely from the loss of angular momentum due to the outflow 
of plasma from a magnetized disk. In this limiting case, the dissipation due to the viscosity and finite 
electrical conductivity of the plasma can be neglected. We have obtained self-consistent, self-similar 
solutions for dissipationless disk accretion. Such accretion may result in the formation of objects 
whose bolometric luminosities are lower than the flux of kinetic energy in the ejected material. 
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Disk-accretion mechanisms can provide clues about a number of processes occurring when matter is accreted onto a 
gravitating center. It is thus no wonder that such mechanisms have been intensely studied for more than 30 years. The 
detection of directed ejections from both the cores of active galaxies [1] and Galactic and intergalactic sources, such as 
Q ' SS433 [2], young stellar objects [3], and Galactic super luminal sources [4], has revealed an unambiguous connection 
between jets and accretion disks. The presence of a directed ejection in a source is essentially always associated with 
evidence for the presence of an accretion disk [3]. One exception is the jet like structures in plerions [5], whose nature 
is obviously different [6]. Thus, the collected observational data indicate that disk accretion is frequently (or even 
always) accompanied by the narrowly coUimated ejection of a substantial fraction of the accreted matter from the 
I— I \ source. This picture gives rise to a number of problems. What is the mechanism that leads to some fraction of the 
■ matter falling on to the gravitating center being ejected from the source? How is this matter collimated? What 
makes the collimated ejection stable enough to propagate to great distances from the source? What accelerates the 
Q , charged particles that generate the observed radiation in the jets? This list of questions is far from complete. Here, 
' we consider just one of them: what mechanism ejects some of the matter falling on to a gravitating center in the 
, course of disk accretion? To fall onto a gravitating center, matter with nonzero angular momentum must somehow 
lose angular momentum. This angular momentum must be carried away by matter, so that disk accretion should 
I inevitably be accompanied by an outflow of matter; the only problem is to determine its rate and direction. In the 
T-H • classical theory of . accretion disks [7], angular momentum of the accreted matter is lost via viscous, hydrodynamical, 
J> ' turbulent stresses. The classical theory does not include a magnetic field. In [8], however, the magnetic field was used 
. to provide high viscosity. This study also considered the possibility of a matter outflow along the axis of rotation of the 
' accretion disk in a hypercritical regime, when the radiation pressure exceeds the gravitation. A quantitative theory of 
~q I such an outflow is given in [9] , which shows that the outflow already begins in a state with subcritical luminosity. This 
• ' ejection mechanism may be important in some sources with hyper-Eddington luminosities; however, it is clear that, 
' in many cases (e.g., in young stellar objects), the ejection is driven by another factor, since the radiation pressure 
in these objects is too small to support this mechanism. Even in SS433, the ejection energy cannot be provided 
purely by radiation pressure [10]. In standard disk-accretion theory, angular momentum is lost via viscous stresses 
due to the differential rotation of the disk. Radiation-efficient disks are geometrically thin, with a local temperature 
that is substantially lower than the virial temperature. It is very difficult to provide an appreciable outflow from the 
disk under these conditions, since the disk matter is in a very deep potential well. Recently, accretion in radiatively 
5^ inefficient disks (ADAFs) for which the ion temperature is close to the virial temperature have been widely discussed 
I [13, 14]. Generally speaking, the matter in such disks can flow out vigorously, with the subsequent formation of a 
jet. However, the existence of ADAF disks is in doubt [15]. Therefore, it remains important to search for alternative 
mechanisms for the ejection of accreted matter. Ejections due to the presence of a magnetic field in the disks are of 
particular interest. 

Blandford and Payne [16] have suggested one of the most promising mechanisms for jet outflows. Matter moves in 
the disk with close to Keplerian velocities. If a magnetic field penetrates the disk, the matter is blown off the disk 
by centrifugal forces if the inclination of the magnetic-field lines to the disk is smaller than 60°. This mechanism is 
particularly easy to understand if we imagine the motion of the plasma along a magnetic line to be like that of a 
bead on a rotating wire. The bead moves in the disk along a circular Keplerian orbit, with the sum of all the forces 
being zero. If a field line is inclined to the equator at an angle smaller than 60°, any shift of the bead along the 
magnetic-field line will result in an increase in the centrifugal force and a decrease in the gravitation force. In this 
case, the position of the bead on the disk is unstable. In addition to forming an ejection, this mechanism can explain 
the loss of angular momentum, which is carried away from the disk by the wind. The influence of the matter outflow 
on the dynamics of the accretion disks themselves was considered in [17-19] in the context of young stellar objects. 
It was shown that the loss of angular momentum due to the wind from the disk can be more efficient than losses due 
to turbulent viscosity. This raises the issue that disk accretion onto a gravitating center should take into account the 
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fact that some of the angular momentum is carried away by the disk wind. Ferrcira [21, 22] attempted to solve this 
problem. His formulation is similar to that used previously by Bisnovatyi-Kogan and Ruzmaikin [20] in one of their 
pioneering studies on the theory of disk accretion in a magnetized plasma. It is assumed that the interstellar plasma 
is the source of an almost uniform magnetic field. In this case, in order for matter to fall onto the gravitating center, 
not only must angular momentum be lost, but matter must also diffuse across the magnetic field. In addition to a 
high turbulent viscosity, this requires a low electrical conductivity of the plasma, which must be associated with a 
high level of turbulence. In other words, dissipation plays a key role in this case. We will use another formulation 
of the problem. Figure la presents the difference between our formulation and that used by Ferreira et al. [21, 22] 
(left). Interstellar matter with some initial angular momentumand a relatively uniform magnetic field is accreted 
onto a gravitating center. The accreted matter diffuses across the field lines, while the forming wind flows out along 
magnetic lines stretched from interstellar matter. Figure lb (right) presents another situation. Matter is accreted onto 
the compact object along with the frozen-in magnetic field, forming the accretion disk. During the accretion, magnetic 
lines emerge from the accretion disk and are stretched. A specific mechanism for this stretching of magnetic lines 
from the disk was proposed and studied numerically in [23]. It results in the formation of a corona and, accordingly, 
a magnetized disk wind. In this case, the accretion disk itself is a source of the magnetic field of the outflowing wind, 
just as the Sun is the source of the magnetic field in the solar wind. Note that precisely such a formulation of the 
problem was used in the study [16] of the centrifugal mechanism for the outflow of matter from an accretion disk. 
Here, we study not only the outflow itself, as was done in [16], but also the influence of the plasma outflow from 
an accretion disk on the accretion. This is not only important for the physics of accretion disks, but is also directly 
applicable to the interpretation of the available observational data. Observations of Galactic superluminal sources 
[4, 24] are of particular interest, since they provide conclusive evidence that plasma ejections substantially affect the 
accretion process. 



II. FORMULATION OF THE PROBLEM 



In general, viscous stresses in accretion disks are undoubtedly important for angular-momentum transport. Here, 
we primarily consider the influence of the wind on the accretion dynamics. Therefore, we will consider the limiting 
case when all of the angular momentum is carried away from the disk by the wind with no dissipation. The possibility 
of such dissipationless accretion can easily be understood based on the following simple consideration. Let one of the 
components of a binary rotate around a fixed gravitating center (Fig. 2). In the absence of dissipation and wind, 
this rotation will continue indefinitely. Let us now suppose that a magnetized wind flows from the rotating object, 
and that the Alfven radius is Ra- We will assume that the wind is spherically symmetrical and its velocity is vq. For 
simplicity, we will use here the mechanical analogy of beads moving along wires for the motion of the plasma in a 
strong magnetic fleld. This kind of motion occurs up to the Alfven radius, beyond which the plasma moves essentially 
freely. We assume that there is no intrinsic rotation of the object, so that its angular momentum is M ■ 1^, where M 
is the mass of the object, and = flB^^ is the specific angular momentum due to its orbital motion. We denote L 
to be the total angular momentum contained in a sphere with radius Ra with its center at the point Ro(i) at time 
t. At time t + St. the angular momentum contained in this sphere (with its center at Ro(t)), will differ from L. The 
equation for the variation of the angular momentum has the form 

^ = - [ra V]^ Vn dS. (1) 

Here p is the density of the plasma, u„ is the velocity normal to the integration surface, the integration is carried out 
over the Alfven surface, and the z axis is directed along the angular-momentum vector. We can see from Fig. 2 that 
Va = Rq + Ra, the velocity of the plasma is v = vqYI + [Jlra], where n = Ra/i?a is the unit vector in the direction 
Ro, and n is the angular-velocity vector for the rotation about the fixed center. After integrating over the directions 
of the vector n, we obtain 

^=-Im^{2RI + RI). (2) 

To interpret this result, note that Lm = —MflR^ is the variation of the angular momentum due to the decrease of 
the mass of the rotating object. It can be shown that the rate of variation of the angular momentum of the wind is 
Ly = —MQRq/3. This can be determined as the difference between the angular momenta Lv{t + 6t) and Lv'(i) of 
the wind particles in a sphere of radius Ra with its center at the fixed point Ro{t). The sum Lm + Ly yields the first 
term in (2). Therefore, the second term describes the rate of variation of the specific angular momentum associated 
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with the orbital motion: 

Thus, the outflowingwind carries away the angular momentum of the rotatingobject. 

This simple example illustrates an important property of accretion in the presence of a magnetized wind. Since 
the accretion is not accompanied by heating, the wind carries away not only angular momentum, but also rotational 
energy. Dissipationless accretion can result in the formation of ejections for a relatively low luminosity of the object. 
This makes this accretion mechanism particularly interesting: it can explain the unusually high efficiency for the 
transformation of the gravitational energy of the accreted matter into the kinetic energy of the jet in the unusual 
source SS433. The bolometric luminosity of SS433 is about 3- 10^^ erg/s [25], while the flux of kinetic energy from the 
object is at least twice this value (> 6 • lO''^ t!rg/s) [12], which is difficult to understand in classical accretion models 
[10]. It is reasonable to consider dissipationless accretion in an ideal magnetohydrodynamical (MHD) approximation. 
This leads to the model in the left panel of Fig. 3: the wind stretches magnetic-field lines from the accretion disk. 
The polarity of the lines varies chaotically, since the total flux of the magnetic field emerging from cither side of the 
accretion disk is, on average, zero. At first glance, it would seem impossible to describe the dynamics of a wind in the 
presence of a magnetic field with arbitrarily varying polarity at the disk surface. However, there is a circumstance 
that simplifies the situation radically. As was noted in [26] , the dynamics of an ideal plasma are invariant with respect 
to the direction of the magnetic-field lines. Therefore, if the direction of the field lines varies so that the polarity 
of the magnetic field is the same on each side of the accretion disk, this does not affect the plasma motion, and 
we are able to consider a self-consistent outflow of plasma from an accretion disk with an azimuthally symmetrical 
magnetic field, shown schematically in Fig. 3 (right). Two features make this approach different from the well-known 
formulation of Blandford and Payne [16]. First, the outflowingwi nd must be consistent with the accretion rate in the 
disk, whereas, in [16], the disk and wind parameters were specified independently. Second, since the plasma moves 
toward the gravitatingc enter, there is a nonzero azimuthal electrical field, ^ 0. Note that a similar formulation 
was given by Contopoulos [27]. However, his study was not continued, since it was clear that, in this case, magnetic 
flux should accumulate at the center. If the polarity of the fleld is the same everywhere, it cannot be annihilated, 
and the accretion will inevitably be halted by the pressure of the magnetic field. As we noted above, this problem is 
removed because we are solving for a plasma flow in the presence of a magnetic field of variable polarity. Therefore, 
the magnetic flux does not accumulate at the center, and is always equal to zero, on average. Since the magnetic-fleld 
lines possess different polarities, they can reconnect, so that the magnetic-field pressure at the center will not increase 
without bound. Finally, we will assume that the plasma flow in the wind is self-similar. This assumption is natural if 
we are interested in the motion of plasma directly above a very thin disk away from its edges. Then, the problem has 
no parameters with the dimension of length, and the flow becomes self-similar. This approximation has been used to 
describe magnetized winds from accretion disks in numerous studies, starting with [16]. Vlahakis and Tsinganos [28] 
present a detailed description and classification of all possible types of selfsimilarity. Note also the study [29] , which 
gives selfsimilar solutions describing nonstationary accretion. 

III. GENERAL RELATIONS 

The dimensionless steady-state equations for an ideal, cool plasma (with pressure P = 0) can bewritten (see, for 
example, [30]) 

p(vV)v = - 1 V + (B V)B - p g , (4) 

div(pv)=0, (5) 
divB = , (6) 

curl[vB]=0. (7) 

Here, the coordinates, density, velocity, and magnetic field are expressed in units of their characteristic values iio, 
POj ^^Oj ^0, respectively; in addition, 

P.vl = f^. (8) 
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The parameter 

g = GM/ (R^vl) , (9) 

where G the gravitational constant and M the mass of the central object. 

Let us consider azimuthally symmetrical solutions for these equations. The dimensioned variables will be chosen 
as follows. Let us arbitrarily select a field line, and adopt Rq as the distance from the center to the point where this 
line crosses the disk. We will take vq to be the Keplerian velocity of a particle at a distance Rq from the center. We 
denote 2 tt to be the magnetic flux through a circle with radius Rq located in the disk and with its center on the 
axis of symmetry. Bq can then be defined by the relation Bq = ^'o/i?o- The dimensioned density Bq = ^q/Rq is 
defined by relation (8); g = 1. Below, we will use only dimensionlcss variables. 

In the azimuthally symmetrical case, the field lines form a two-parameter family of curves, which can be presented 
in the form 

'R{a,ijj,ip) = {r{a,ip) cos r{a,ijj) sin (j), z{a,ip)) , (10) 

where 4> = ip + ri{a, tp). The parameters ^, if specify a field line, and a varies along a given line (B ~ dR/da). The 
substitution (p ^ (p — (po is equivalent to a rotation by (p>o of the field line as a whole. The functions of two variables 
r{a,'ip), z{a,tp), and ri{a,'ip) are to be determined from the solution of the MHD equations. Following [31, 32], we 
will call the parameters a, tp, and (p the frozen-in coordinates. The cylindrical coordinates (r, (J), z) can be expressed 
in terms of the frozen-in coordinates by the formulas 

r = r(Q!,V'), (j) = ip + 'n{a,ip), z = z{a,ip). 

Let us make a transformation from Cartesian to frozen-in coordinates in (4) (7). J denotes the Jacobian of the 
transformation: 



d{X. Y, Z) 
d{a,tp,(p) 



where the subscripts denote differentiation with respect to the corresponding variables (R^ = dTl/da, . . .). Since 
field lines do not cross, there is a one-to-one relation between the frozen-in and Cartesian coordinates, so that the 
Jacobian J never vanishes. Note the following useful equalities, which can be used to easily change between Cartesian 
and frozen-in coordinates: 

R„V = A, R,V = A, R,V = A. (12) 



and also 



and 



^ = [V^,V^], ^ = [V<^,Va], ^ = [Va,V^], (13) 



l3tlM=ya, l?^^=V,p, [5^!iM = Vv. (14) 



The operator V in frozen-in coordinates has the form 



V = - I [R^R^] ^ + [R^R J ^ + [R„R^] ^ j . (15) 
Suppose that the magnetic field strength is specified by the equality 

B = ^ = [VV,V<p]. (16) 

In this case, (6) becomes an identity (it essentially expresses the field in terms of the Euler potentials). Expression (16) 
implies that we take the magnetic flux (in units of 27r^'o) through a circle with radius r as the frozen-in coordinate 
Note that the ratio R^/J is invariant with respect to the reparameterization a — > a' = g{a,tl),(p); i.e.. 
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where g{a,ip,ip) is an arbitrary function and J' = d{X,Y, Z)/d{a' ^tp-ip) = J/ga- This can be used to assign an 
arbitrary value to the Jacobian J without varyingthe field strength. It is convenient to relate J with the plasma 
density: 

J=l/p, (17) 

in which case the magnetic field is 

B = pR«. (18) 
Let us introduce the three mutually orthogonal unit vectors 

e^'^) = (cos 0, sin 0,0), e^"^) = (- sin0, cos </>, 0), 6^ = (0,0,1). (19) 

We obtain from (10) 

R„ = Ta e^"") + r 77a e^"^) + z„ e(^) , \ 

= rv,eM+rry^e(<^) + ^^eW, \ (20) 
R^ = re(^). J 

Therefore, the Jacobian of the transformation is 

J = r {zocT^^, - z^ra) , (21) 

and, in accordance with (17), the plasma density is 

P = ^ r • (22) 

Let us expand the plasma velocity into three noncoplanar vectors: 

V = /Ra + ORy + £R^ = 

= ifra + er^) e^^^ + {fz„ + ez^) e^^) + r (/r?« + er]^ + Q) e^'^) , (23) 

where the functions /, fJ, e may depend on the frozcn-in coordinates. It turns out that, in both the stationary and 
nonstationary cases [31, 32], Eqs. (5) and (7) in frozen- in coordinates can be integrated in general form. Substituting 
(23) and (17) into the continuity equation (5) and transforming to frozen-in coordinates, we obtain 

df on de _^ 
da dip dip 



The electrical-field strength is 



E = --[vBl = -(eVy-OVV'). (25) 

c c 



It follows that (7) in frozen-in coordinates has the form 

curl[v B] = [VO W] - [V£ V<^] = p (^1^ R^ + 1^ Rv, - + ^) R.) = . (26) 

The derivative d^l/dip due to the assumed azimuthal symmetry. Solving(24) and (26), we find that e = const, while 
the functions / and may depend only on the single variable ip. In this case, (5) and (7) are identically satisfied. The 
gradientVV' has only r and z components; the first term in (25) specifies the azimuthal component of the electrical 
field: 

E^ = —. (27) 
cr 

The velocity of the accretion disk u is also expressed in frozen-in coordinates. Let us suppose that the disk surface 
corresponds to the parameter a = 0; i.e., z{0,ip) = 0. Then, the derivative Zjp{0,'tp) = 0. The plasma density above 
the disk surface is 

Po = ^— , (28) 

rz^r^ 



The z component of the magnetic field is 



B.o = P^a|o--^ . (29) 



The radial component of the electrical field at the disk surface is 

^rO = --[vB]X= — L(n + e^)| (30) 

while the azimuthal component is defined by (27). Using the frozen-in condition 

E = -i[uB], 

the continuity of the z component of the magnetic field at the boundary between the plasma and disk, and the 
continuity of the r and (l> components of the electrical field, we obtain the components of the disk velocity 

Ur=er^|g, Zi^ = r (il + er;^)!^. (31) 

After transforming to frozen-in coordinates, the Euler equation (4) assumes the form 

= -\ {[R*R-v']^ + [R^Ra]|^ + [RaRv]^|(pRa)'-9;^- (32) 

Substituting R the form (10) and p in the form (22), we can obtain equations for the functions r(a,z/'), z(a, V'), and 
r/(Q!,V)- 

IV. BOUNDARY CONDITIONS 

Let us study the relation between the field components and the parameters of the accretion disk, assuming the disk 
is infinitely thin (which is valid for a cool plasma). We will specify the volume density of the disk in the form 

Pi{r,z)=a{r)b{z), (33) 

where 6{z) is a delta function and a(r) is the mass surface density. Equation (4) is valid for ^; > (or < 0); i.e., 
beyond the disk. The dynamics of the plasma in the disk are specified by the equation 

Pd(uV)u + p(vV)v = -i V B2 + (BV)B - (p, + p) g . (34) 

Wc will integrate (34) over z in the small interval (— zo, zq) and take the limit z^^ — > +0. This operation will be 
applied to each term of the equation. In cylindrical coordinates, 

V = e«|-+e(^)|-+e(-)47, (35) 
or oz rocp 



so that 



Pd (uV)u = a{r) 5{z) [^r^ + u^-^^+u^^] (e^'^^Ur + e^^^u^ + e^"^)^^) = 
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When deriving(36), we took into account azimuthal symmetry, the equahty Uz = 0, and the relations de'^'^^ /d<f> = e^"^\ 
de^'^'>/d(j) — —e^^\ Integrating (36), we obtain 

JpAnV)udz^a{r) - ^J^^j + + ^) e(*) . (37) 

-Zo L \ / J 

If the magnetic field and the plasma velocity above the accretion disk (z i 0) are known, then the solution in the 
region under the disk continues in accordance with the rules 

Bz ^ Bz 1 Br —Br , B^ —B^ . (38) 

Vz -Vz , Vr ^ Vr, ^ . (39) 

Since the field components B^, B^ at the equator (for z = +0), this implies that the field line is discontinuous there. 
Using relations (39), it can be shown that, in the 



J p{yV)\dz = 0. 



Next, we obtain 

j VB^ dz^Q, (40) 

-Zo 

since B^ is a continuous function of z. 

The term (BV)B will be written in the form 

(BV)B = + Bz^ + B,^^ (i3„e(^^) + S.e(^) + B.e^^)) . (41) 

If the integration is taken over an infinitely small interval, a nonzero contribution will be obtained only from terms 
containingth e derivative d/dz (a derivative of a discontinuous function yields a 5 function). After the integration, 
we obtain 

Zo 

J{B\/)Bdz = 2Bz{Brer + B^e^). (42) 

-Zo 

The integration of the last term in (34) is trivial due to the presence of the 6 function: 

Zo 
-20 

Collecting these various relations, we obtain two equalities (for the r and 4> components in the equation), which are 
a direct consequence of (34) in the model with an infinitely thin disk: (|34p : 

dur 2 q 

- — = - BrBz - 4 , 44 
or r a 

-^{tu^)^-B^Bz. (45) 
r or a 

The conservation of mass yields another relation: 

--^iraur) + 2jz^0, (46) 
r or 

where jz — pvz is the density of the z component of the mass flux. In (44) (46), the magnetic- field components 
and plasma parameters are taken for the disk surface {z = +0). Equations (44)-(46) play the role of boundary 
conditions, which must be satisfied by the solution. They were derived from the equations of motion in differential 
form. Obviously, the same conditions can be obtained using the integral form of the equations of motion. 



8 



V. SELF-SIMILAR SOLUTIONS 



We are primarily concerned here with the plasma dynamics immediately above the disk, at distances z much smaller 
than the radius of the disk i?disk- In the limit z <§; i?disk: only two parameters with the dimensions of length remain 
in the problem; z and r. Therefore, the solution will be self-similar in this limit [11]. Let us underscore an important 
feature of these solutions: they describe flows only at small distances from the disk, and are not applicable at distances 
comparable to or exceeding the size of the disk. Equations (4)- (7) have self-similar solutions of the form . 

v(r, z,(^) = r"''" v(z/r. (i) . ^ 

p{r,z) = r-^opiz/r), \ (47) 
'B{r,z,(l)) = r-^''B{z/r,(l)) . J 

The superscripts (5„, Sp, and 6b are determined from the following conditions. Substituting (47) into (4) leads to 
the equations 



Equality (27) is satisfied under the condition 
Thus, 

and hence, 



26b -6p = 26y = l. (48) 
6^ + 6b = 1. (49) 



(5„ = (5b = - , 6p = (50) 



v(r, z,0) = r ^/^v(z/r, 0) 

p{r,z) = p{z/r), ) (51) 

B(r,^,0) = r-V2B(0/r,0). 

Note that it is possible to determine the subscripts unambiguously only for Ur ^ 0; i.e., when there is advection 
of the plasma toward the center of gravitation. When Ur = 0, one of the subscripts is arbitrary, so that there is a 
single-parameter family of selfsimilar solutions. A large number of such solutions have been studied previously (see, 
for example, [28] and references therein). For clarity, we will compare our results with those of [16], which is similar 
to our study in terms of the formulation for the wind. The case 6p = 3/2, 6b = 5/4. was considered in [16]. 
In frozen-in coordinates, the law of similarity (51) is described by the relation 

R(a,V,^)=V''/'R-(^,<p) , (52) 

Therefore, the two-variable functions in (10) are expressed in terms of one-variable functions. Denoting s = a/'ij}, we 
can write 

r{a,i)) = ijj'^/^ f{s) , z{a,il}) = ip"^/^ z{s) , ri(a,ilj) = fj{s) . (53) 
Let us verify the equivalence of (51) and (53). Substituting(53) into (22), we derive the plasma density: 

p{s) = I ^ , (54) 

where the subscript s denotes differentiation with respect to s. This function depends only on s, or, if we switch to 
cylindrical coordinates, on the ratio z/r, as in (51). The components of the magnetic-field strength are 

B^ = ^-V3pf,, B^=V-'/3p^., B^ = ^-^/^pffi,. (55) 
We find for the components of the plasma velocity 

Vr = ^IJ-^/^(^^sf + {f-es)ny (56) 

V. = iJ-'/^(^^sz + {f-es)z,y (57) 

= i^-^'^f{^n+{f-es)fis)- (58) 
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The law of similarity v r will be satisfied if / = const, while O depends on tp as follows: 

= (l/ip . (59) 

represents the angular velocity of the selected field line in units of the Keplerian angular velocity. Substituting(52) into 

(32). wc obtain a system of three ordinary second-order differential cquations"'^ for the functions r(s), 5(s), and fi{s). 
Since second derivatives appear in the equations in linear form, the system can be solved for the second derivatives. 
For £ 7^ 0, the right-hand sides of the equations depend explicitly on s; f and s appear in the equations only in the 
combination / — es. The parameter s has an arbitrary zero point. Indeed, shifting the reference point for s via the 
substitution s — s' + sq in equivalent to redefining the constant /: f ^ f = f — ssq- Wc will take the surface of the 
accretion disk as the reference point for s: z{s = 0) = 0. 

The resulting equations can be written as a system of six ordinary first-order differential equations. We introduce 
the six-dimensional vector with the components 

^l=fs, £.2 = Zs, 6 = % , (60) 
?4 = r, £,5 = Z, £,6 = fl. (61) 



Then, for i = 1, 2, 3 the equations have the structure 

^^1,2,3 _ -^1,2,3 

ds ~ D 

while for i = 4, 5, 6 



(62) 



ds 



= 6,2,3 . (63) 



The denominator D is the determinant of a certain 3x3 matrix; it appears when the initial system of equations is 

solved relative to the highest derivatives. The functions Ni and D depend on ^i, . . . , ^5; ^6 docs not appear in the 
equations due to the azimuthal symmetry of the problem. In addition, the right-hand sides of (62) depend on s . For 
fixed ij) and if, the curve R = tj)^/^ R(s, ip) specifies a magnetic-field line in parametric form. In 

the case of self-similar solutions, it will suffice to find one field line; others can be derived using the similarity 
transformation (multiplication by ^p^^^) and rotation about the z axis by some angle (p. The field line corresponding 
to the frozen-in coordinate ip = 1 crosses the disk at the distance r(0) = 1. The value 7y(0) can be arbitrary, for 
example, zero. The solution of the system of differential equations should then satisfy the initial conditions 

f(0) = 1 , z{0) = , ?7(0) = . (64) 

The frozen-in coordinate tp and cylindrical coordinate r are related by the simple expression 

r = ^^/^ . (65) 

To obtain an unambiguous solution of the system of equations, we must have another three conditions. 

It follows from (54) that the plasma density at the disk surface is p = 3/(2 5^°^) Therefore, the components of the 
magnetic field in the disk are 

q ~(0) q q ~(0) 

^^°^ = ^5o)V'-^'^ 4°^ = ^V'-^/^ < = ^^^-^/^ (66) 

Zs Zs 

where the superscript denotes values for z = -1-0. 

Let us write the boundary conditions (44) (46) for the self-similar solution. Substituting(65) into (31) and applying 
the initial conditions, we obtain the velocity of the accretion disk: 

W^ = ^eV'"'/^ W^ = fiV-^/^ (67) 



The equations are cumbersome, and we do not present them here. The analytical manipulations were performed using the Maple 
package. 
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The velocity component Ur < 0, so that s should also be negative. Together with the plasma velocity. and 
depend on r as r~^/-^. It follows from (44) or (45) that the mass surface density of the disk is a linear function of r: 

a = fjbr , (68) 

where = const. In (46), jz is expressed using (54) and (56) as follows: 

j. = \fr"'. (69) 

We derive from the continuity equation (46) 

^l = Z^/\e\. (70) 
Substituting(66), (67), and (68) into the boundary conditions (44) and (45), we obtain 

«?*=-5/!j4°'- (72) 
This last expression can be used to calculate the azimuthal component of the magnetic field at the equator: 

Bf = -\fnr"'. (73) 

We can see from these conditions that the rotational velocity of a field line is lower than the Keplerian velocity. Only 
in the limit £ ^ do we obtain f2 = 1, which corresponds to a Keplerian accretion disk. Note that the boundary 
conditions in the disk fully specify the azimuthal magnetic field above the disk. Since e does not appear in (72) and 
(73), these equalities should be satisfied even in the limit e ^ 0. 

The conditions (44) and (46) make it possible to relate in a self-consistent way all the parameters of the accretion 
disk and of the outflowing plasma, and, in particular, to take into account the decrease in the surface density of the 
disk due to the expansion of the plasma. The equality (45) leads to the additional condition (72), which must be 
satisfied by the solution. 



VI. SOLUTION IN THE VICINITY OF THE DISK 

Let us consider an important limiting case, when the radial velocity of the disk is small compared to the Keplerian 
velocity, |£| <C 1. Let us also assume that the velocity of the plasma outflow from the disk is fairly low. Then, the 
following equations are valid in the vicinity of the disk: 

— =rsG, —=ZsG, — =r)sG. (74) 
as as as 



Here, 



6/3 5^(f2+52+^2) ' I ) 



where 5f = f — 1. To derive (74), Q in the form (71) is substituted into the exact equations (62) and expanded in 
small parameters. Equations (74) are valid when 5f and z are small compared to unity and 

«1, /^z, «1, /^r?, «1. (76) 

These conditions follow from the requirement that subsequent terms of the expansion be small. It follows from (74) 
that 



:^(o)exp(yGds), (77) 
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where is any one of the three functions fs, Zg, or fjs. After integrating (77) with the initial conditions ^(0) = 0, 
(5f(0) = 0, 77(0) = 0, we find that a field line is straight in the vicinity of the disk: 

(5f(s) = kr z{s) , 77(5) = kn z{s) , (78) 

where kr and fc,, are constant. As we can see from (72), fc^ = —2 //9 in the approximation considered, kr specifics the 
direction of the poloidal field in the vicinity of the disk. Using equality (78), the second of Eqs. (74) can be presented 
in the form 



where 



"^'^ =w'z-\e\C, (79) 



1 9fc,^-2/fc, 

^ 6,p A;2 + A-2 + l • ' 

When w"^ < 0, (79) has an oscillating solution, which is physically unacceptable. The inequality k^ > 1/3, obtained 
in [16], must be satisfied by field lines moving away from the disk. This implies that a magnetic-field line must be 
inclined to the disk at an angle smaller than 60°. 

Unlike [16], we also obtained another constraint. When > 0, the solution (79) has the form 

lei C z^°^ 
= (1 - cosh(u;s)) + — sinh(«;s) , (82) 

where Zs"-* is an arbitrary constant equal to the derivative Zs in the disk. It follows from (82) that solutions cannot 
exist for arbitrarily small zi*'' . The field line will move away from the disk only when 

5(0) 5 = Mi 



>Zsmin=^. (83) 



This is equivalent to a constraint on the z component of the plasma outflow velocity: 

Vz\^_0>V.min = f—- (84) 

The nature of this constraint, which appears when e ^ 0, is easy to understand. In this case, each particle of the 
plasma moves along its orbit with a velocity that is smaller than the Keplerian velocity. In the z direction, the particle 
is in the effective potential 

U{z) = ~w'^ z^ + \s\Cz. 

In order to be ejected from the disk, the particle must overcome the potential barrier corresponding to the difference 

between the gravitational and centrifugal forces. 
The plasma density in the vicinity of the disk is 

"'^.'TU^iZZ^ ^' '''' 

In the disk itself, it is finite: 



2 I Zs ch( ws) sh(t(;s) 

w 



I 3 3w 



In the limit £ ^ 0, zi^^ 0, the density in the disk increases without bound, as in [16]. 
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VII. PASSAGE THROUGH CRITICAL SURFACES 

The general solutions of the equations describing stationary MHD flows are singular on certain surfaces, usually 
called critical surfaces. These have already been found for onc-dinicnsional flows, as was proposed by Parker for the 
solar wind [33]. In one dimensional solutions, they are manifest as critical points, while, in the three-dimensional 
case, a set of critical points forms a surface. In the one-dimensional solutions, the critical points coincide with sonic 
points, where the local velocity of weak perturbations becomes comparable to the flow velocity, while the type of 
equation changes from elliptical to hyperbolic and back. Even the earliest solutions obtained for axisymmetrical flows 
indicated that this is not the case in general [16]. The nature of critical sm-faccs and their role in the formation of 
stationary MHD solutions was investigated in [34, 35]. It turns out that critical surfaces divide regions of solutions 
with different causeeffect relations, and these surfaces do not, in general, coincide with the surfaces where the form 
of the equation changes. 

Critical surfaces as applied to self-similar solutions were studied in [36]. The parameters of critical surfaces for 
self-similar flows with ^ Q were studied in detail in [27]. In these solutions, the locations of the critical surfaces 
arc specified by the zeros of the denominator D in (62). In the general case (with a nonzero wind temperature), 
the denominator vanishes at the slow magnetoacoustic separatrice surface, Alfven surface, and fast magnetoacoustic 
surface. Regularizing the solution on these surfaces makes it possible to determine the density of the plasma flux 
from the disk and the two tangential components of the magnetic field at the disk surface (the normal component is 
specified by the conditions adopted for the problem). In the case of a cool plasma, the slow magnetoacoustic velocity 
is zero, and one of the regularization conditions disappears. Therefore, the plasma flux from the disk surfacemust be 
specified for a cool plasma. This means that the ratio f of the density of the plasma flux from the disk and the z 
component of the magnetic field is also specified. A similar approach was used in [16] (where k is analogous to /). 

Thus, in self-similar flows of cool plasma, regularizing the solution on the Alfven and fast magnetoacoustic crit- 
ical surfaces specifles both the azimuthal and tangential components of the magnetic field based on the fiow. One 
characteristic feature of self-similar solutions is that, as a rule, it is not possible to regularize them on the fast magne- 
toacoustic critical siirface, only at the Alfven surface. Recently, Vlahakis et al. [37] were able to regularize a selfsimilar 
solution on the fast magnetoacoustic surface, but, even in this case, the solution disappears very rapidly behind this 
surface. Therefore, as in [16], we will not try to regularize the solution on the fast magnetoacoustic surface. As 
a result, one of the components of the magnetic field (the azimuthal component in [16]) remains free and must be 
specified a priori. The other component is determined from the solution and the regularization condition at the Alfven 
surface. The denominator D in (62) contains the factor 

e^=p-{f-ssf. (87) 

In the disk, ei > when pz'f'^ < 3/2. Since the density decreases along a field line, ei vanishes at some point on 
the Alfven surface. To keep the solution finite, the functions Ni must also vanish at this point. After substituting 
p = (/ — £s)^ into Ni, all the functions Ni possess the common factor 



+ I (s f^ze^fjs + 8 fz^e^rj^ + 12 fz^s Cl 2^-12 f^e ^z^ + ^Q f^zn^fjs^ - 108 j fs+ 



+ {-3Qf%h^f)s -Sr^z^i^fis - 12 r^zeClz"^ - Sr^z^e^risz'^ - 27 fm - 27 fj.zfj 

+ 54 - 108 bn fm , (88) 

where & = / — es. Therefore, the condition that the denominator and three functions Ni in (62) simultaneously vanish 
results in the system of two equations 

ei = , 62 = , (89) 

which can be solved analytically for fs and Zs- Consequently, the derivatives [fs, ^s) can be taken to be specified at 
the critical point (f, z). 
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VIII. SOLUTION IN THE LIMIT e 0. 



It is of interest to consider the solution in the Umit of very slow accretion, when e — > 0. This case is of interest 
because the disk becomes Keplerian in this limit (17 = 1), so that our results can easily be compared to those of [16], 
and the physics of accretion made possible by a wind outflow becomes particularly simple. 

In the case of small e, es can be neglected compared to /. The parameter s then disappears from the right-hand 
sides of (62), and the simplified equations have two integrals of motion: the momentum L and the energy W. For the 
self-similar solution, 

L = f^(^fn-fjs{p-f)) , (90) 



w 



l /■2 /~2 I -2-2 , -2\ l /S2-2 ^ 

-/ (r„ +r 71^ + z.) — -it r , 

2-' V « ' Is < sj 2 ^y~T^rf2 



(91) 



where we can assume f2 = 1 on the right-hand sides of the equality. 

Substituting the values at the disk surface into (90) and using (64) and (72), we find 



(92) 



Since zi°^ > 0, L < 3 /. This new important constraint on the angular- momentum flux carried away from 

the disk follows from the boundary conditions at the disk. It means that a relatively small flux of angular momentum 

must be carried away by the wind to enable accretion. In addition, the boundary conditions at the disk specify the 
azimuthal component of the magnetic field, and the problem is fully specified even when the solution is regularized 
at one Alfven point. 

As we can see from (84), in the limit e 0, the initial velocity of the plasma outfiow can be arbitrarily small. In 
this case, the energy tends to the minimum W^ain — —3/2 [see (91)]. Assuming that the conditions (76) are satisfied 
at the disk surface, we obtain L = 4/3 f,W = —3/2. Therefore, at the Alfven point. 



„ L 4 
r^ = . = 3. 



/ 



(93) 



The energy integral and Eqs. (89) can be used to express analytically the derivatives fg, Zs, and fjs at the Alfven point 
in terms of z. A straightforward but fairly cumbersome analysis shows that the solutions (89) exist for L = 4:/3f, 
W = -3/2 if 



/> ^^6(5 + 3^3) « 4.4. 



(94) 



z must be confined within the interval 



< Z < Zr 



where 



216/v'16/2- 45 - 736/2 _ 1215' 
1200/2 



1/2 



when / < 8.27, and 



2 ^2351/4 - 159894/2 - 59049 
71 329/2 + 243 



(95) 



(96) 



(97) 



when / > 8.27. It follows from (96) that z < 0.327 at the Alfven point; this means that the Alfven point cannot be 
far from the disk surface. 

We solved (62) and (63) in the limiting case e — > numerically using standard techniques. We started from the 

Alfven point, with / as the input parameter. As the initial condition, we took the Alfven point with f = 2/\/3 
and z from the interval (95); the derivatives were selected using the described procedure. At the Alfven point, the 
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right-hand sides of (62) were determined using IHopitales rule. The initial value for z was selected so that the field 
line crossed the disk at the point f = 1. We adopted a value slightly larger than —3/2 for W . The reason for this is 
that, when W = —3/2, the plasma velocity at the disk surface vanishes and the plasma density at the disk displays a 
nonintcgrable singularity: p = 3/(2 wz). Therefore, wc chose W so that .z^f"^ ^ 10"'^ in the disk, or, in other words, 
so that the velocity of the plasma outflow relative to the disk was of the order of / • 10~^ of the Keplerian velocity. 
The requirement that zi^^ be nonzero also results from the fact that the inequality (84) should be rigorous. 

Figure 4 presents the / dependence of z at the Alfven point together with ^min and -J max- The minimum / for 
which the solution exists is /min = 5.28. Zmax = za when / = /min- 

Figure 5 presents a line of the poloidal magnetic field for / = 6, which does not differ much from /min- kr = 0.96, 
so that the line is inclined to the equator at an angle that is close to 7r/4; at the disk surface, Zg = 10~^. There is 
no significant acceleration of the plasma; at the Alfven point, Zg — 0.033 and fs — 0.010, so that the components of 
the plasma velocity Vz and Vr are small compared to the Keplerian velocity. Figure 6 presents the plasma density 
on a field line for various heights. In the case of motion from the disk toward the Alfven point, s varies from to 
Smax ~ 30. Therefore, es can be neglected compared to / when 

|e|«/ASmax. (98) 

This condition, as well as the inequality (76), ensures the applicability of the approximation used in the case considered. 
The analytical solution (82) coincides with the numerical solution along nearly the entire field line, with the exception 
of the immediate vicinity of the Alfven point. 

IX. DISCUSSION 

The conclusion that the plasma outflow from an accretion disk can exert a substantial effect on the dynamics of 
the accretion disk is not new [19]. The originality of our work is our consideration of the case when plasma falls onto 
the gravitating center along with the frozen-in magnetic field, rather than seeping across magnetic-field lines (see, for 
example, [21] and references therein). Our results indicate that dissipationless disk accretion of the type considered is, 
indeed, possible, as is demonstrated by the derived self-consistent solution. The mechanism for carrying away angular 
momentum that we have considered should be taken into account in studies of accretion in sources displaying violent 
ejections of matter. This is particularly important for our understanding of the processes occurring in a number of 
peculiar sources, prominently represented by SS433 [2]. Objects in which disk accretion occurs primarily due to the 
carrying away of angular momentum by an outflowing wind can have bolometric luminosities that are appreciably 
lower than the kinetic-energy flux in the outflowing wind. This is what leads to the most important difflculties in 
explaining the plasma-ejection mechanism in both SS433 and young stellar objects. The derived self-consistent solution 
possesses a number of interesting properties. The boundary conditions (71), (72) fully specify the azimuthal magnetic 
field in the disk. In this case, with fixed / and L, the toroidal and tangential components of the magnetic field are 
determined unambiguously, unlike in the study of Blandford and Payne [16]. The flow turns out to be fully specified. 
This begs the question of whether there exists the fundamental possibility of passing through the fast magnetoacoustic 
critical point, since, at first glance, no free parameters remain. In fact, there is still one free parameter: e, which 
specifies the radial velocity of the plasma in the disk. This parameter can be used to regularize the solution at the 
fast magnetoacoustic critical point. If this is possible, the entire solution will prove to be unambiguously specified. 
An unexpected feature of the obtained solutions is that, in the limiting case of dissipationless accretion, 100% of the 
infalling matter is ejected back out from the disk. In this respect, our results are similar to those of Lery et al. [38], 
who also considered dissipationless accretion, but without the formation of a disk. Both of these results indicate 
that, in any case, dissipation is probably unavoidable, to provide for the infall of some fraction of the matter onto the 
gravitating center. On the other hand, this feature of the flow may well be due to its self-similarity. Additional studies 
arc required before wc will be able to draw final conclusions. A natural extension of our study is the consideration 
of accretion due simultaneously to matter outflow from the disk and dissipation. Even at this stage, it should be 
possible to determine the relationship between the disk luminosity and the flux of outflowing matter, which can be 
compared to observations. 
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Figure 1: Role of the magnetic field in different formulations of the problem. Left: plasma diffuses across the magnetic-field 
lines. Right: plasma falls onto the gravitating center along the magnetic field. The wind flowing from the disk extends some 
field lines to infinity and carries away some fraction of the angular momentum of the infalling material. 




Figure 2: Loss of angular momentum via a magnetized wind. The gravitating center is at the point A. 




Figure 3: Independence of the plasma dynamics of the direction of the magnetic-field lines in an ideal MHD approximation. 
The left panel shows a realistic structure for the flux of the poloidal magnetic field; the total fiux of the poloidal magnetic field 
is zero. In the right panel, the direction of the field lines has changed, while the plasma dynamics remain unaltered. 
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Figure 4: Values za of z{f) for which a field line emerging from the Alfven point crosses the equator when f = 1. Zmax and 
Zmin, which appear in (95), are also indicated. 
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Figure 5: A line of the poloidal field for / = 6. The circle marks the position of the Alfven singularity. 
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Figure 6: z dependence of the plasma density on the field line shown in Fig. 5. 




